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Abstract 

The first observation of XcJ (J=0)l)2) decays to AA is reported using ijj{2S) data collected with 
the BESII detector at the BEPC. The branching ratios are determined to be B{xcO = 
{4:.7t\i ± 1-0) X 10-^ , Bixci AA) = (2.61^;!] ± 0.6) x 10"^ and B{xc2 AA) = (3.311 J ±0.7) x 
10"*^. Results are compared with model predictions. 

PACS numbers: 13.25.Gv, 14.40.Gx, 12.38.Qk 
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I. INTRODUCTION 



It has been shown both in theoretical calculations and experimental measurements that 
the lowest Fock state expansion (color singlet mechanism, CSM) of charmonium states is 
insufficient to describe P-wave quarkonium decays. Instead, the next higher Fock state 
(color octet mechanism, COM) plays an important role [1, 2]. Our earlier measurement [2] 
of the total width of the XcO agrees rather well with the COM expectation. The calculation 
of the partial width of XcJ PP, by taking into account the COM of XcJ decays and using 
a carefully constructed nucleon wave function [3], obtains results in reasonable agreement 
with measurements [4]. The nucleon wave function was then generalized to other baryons, 
and the partial widths of many other baryon anti-baryon pairs predicted. Among these 
predictions, the partial width of Xcj AA is about half of that of Xcj PP (J=l,2) [3]. 

In this paper, wc report on an analysis of the '-fir^TT'pp final state produced in ip{2S) 
decays. Evidence for the decays of Xcj to AA is observed for the first time. The data 
used for this analysis were taken with the Beijing Spectrometer detector (BESII) at the 
Beijing Electron Positron Collider (BEPC) at a center-of-mass (CM) energy corresponding 
to M0(25). The data sample corresponds to a total of about 15 miUion ■0(2-5') decays. 

BES is a conventional solenoidal magnet detector that is described in detail in Ref. [5]; 
BESll is the upgraded version of the BES detector [6]. A 12-layer vertex chamber (VTC) 
surrounding the beam pipe provides trigger information. A forty-layer main drift chamber 
(MDC), located radially outside the VTC, provides trajectory and energy loss (dE/dx) 
information for charged tracks over 85% of the total solid angle with a momentum resolution 
of ap/p — 0.0178-\/l +p'^ {p in GeV/c) and a dE/dx resolution for hadron tracks of ~ 8%. 
An array of 48 scintillation counters surrounding the MDC measures the time-of-flight (TOF) 
of charged tracks with a resolution of ~ 200 ps for hadrons. Radially outside the TOF system 
is a 12 radiation length, lead-gas barrel shower counter (BSC). This measures the energies 
of electrons and photons over ~ 80% of the total solid angle with an energy resolution 
of ue/E — 21%/\/E {E in GeV). Outside the solenoidal coil, which provides a 0.4 Tesla 
magnetic field over the tracking volume, is an iron fiux return that is instrumented with 
three double layers of counters that identify muons of momentum greater than 0.5 GeV/c. 

A Monte Carlo simulation is used for the determination of the mass resolution and de- 
tection efficiency, as well as the estimation of the background. For the signal channels. 
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'il^{2S) — > 7XcJ, XcJ — AA, the angular distribution of the photon emitted in the ■0(25') 
decay is assumed to be that for a pure El transition. The A in the XcJ CM system and the 
daughter particles in the A CM system are generated isotropically. A total of 10000 events 

are generated for each Xcj state with A Ti^p and A vr+p. For the estimation of the 
number of ip{2S) events and the estimation of the systematic error, ip(2S) — >• tt+tt" J/'0, 
J/ip pp events are generated, where the tt+tt" invariant mass is distributed as measured 
in Ref. [7]. 

The simulation of the detector response, including interactions of secondary particles 
in the detector material, uses a Geant3 based package SIMBES. Reasonable agreement 
between data and Monte Carlo simulation is observed in testing various channels, including 
e+e" —>■ 76+6" (Bhabha), e+e~ — >• J/ip — >■ pp, and ip{2S) — >• TT~^7T~J/ip, J/ip —>■ . 



II. EVENT SELECTION 

The analysis uses the same photon selection and charged particle identification (ID) 
criteria as were used in Ref. [8]. When selecting photons it is necessary to remove photons 
produced by hadronic interactions of charged tracks with the detector material. This is 
achieved by cutting on the angle between the neutral cluster and the charged track in the 
BSC. The number of photon candidates in an event is not limited. 

Both TOP and dE/dx information are used for charged particle identification. Proba- 
bilities of a track being a pion (ProbT^), kaon (ProbK), or proton (Probp) are assigned to 
each charged track. Por the decay channel of interest, the candidate events are required to 
satisfy the following selection criteria: 

1. Each charged track is required to be well fit to a three-dimensional helix and be in the 
polar angle region | cos^mdcI < 0.8. 

2. The number of charged tracks is four with net charge zero. 

3. The two lower momentum positive and negative charged tracks are assumed to be the 
TT"^ and the 7r~, and the other two tracks are regarded as the proton and the antiproton. 
Pour-constraint kinematic fits to the decay hypothesis are performed with each of the 
photon candidates, and the one with the smallest is taken as the real photon. The 

probability of the fit is required to be greater than 1%. 
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4. The particle identification assignment of each charged track must satisfy ProhT^ (for 
7r='=) or Probp (for p or p) > 0.01. 

A four-constraint fit assuming 'i/j{2S) —>■ n^n'pp is also performed to select '4>{'2S) AA 
and ■0(2-5') — > t^'^t^~ J /'^, J /'^ — > pp events for checking the reliability of the analysis of 
Xcj — ^ AA and to calculate the total number of t/^i^S) events. The selection criteria used 
are the same as for XcJ A A except that no photon information is used. 

III. EVENT ANALYSIS 

Fig. 1 shows a scatter plot of the Tr+p versus the tt~p invariant mass for events with 
n'^n'pp mass between 3.38 GeV/c^ and 3.60 GeV/c^. The cluster of events in the lower left 
corner shows a clear AA signal. 
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FIG. 1: Scatter plot of 7r~^p versus 7r~p invariant mass for selected 'yTr~^7r~pp events with the 
Tr~^Tr~pp mass in the XcJ mass region. 



Selecting events in XcJ mass region and requiring the mass of 7r~^p (7r~p) to be smaller 
than 1.15 GeV/c^, the n'p (n'^p) mass distribution shown in Fig. 2 is obtained. A clear 
A signal can be seen, and the background below the peak is very small. A fit of the mass 
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distribution gives mA = (1114.6 ± 0.6) MeV/c^, in agreement with the world average [4], 
and a mass resolution of (6.3 ± 0.6) MeV/c^. 
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FIG. 2: Mass distribution of n^p (vr^p) recoiling against a A (A) (mass < 1.15 GeV) for events 
in the Xcj niass region. Dots with error bars are data and the histogram is the Monte Carlo 
simulation, normalized to the A signal region (two entries per event). 

After requiring that both the vr+p and the 7i~p mass lie within twice the mass resolution 
around the nominal A mass, the AA invariant mass distribution shown in Fig. 3 is obtained. 
There are clear Xco, Xa, and Xc2 — * A A signals with low background, estimated using A 
mass side band events. The highest peak around the ip{2S) mass is due to iIj{2S) AA 
with a fake photon. 

Fig. 4 shows the energy deposited in the BSC of the proton track versus the antiproton 
track for events selected as Xcj AA. Since the antiproton will frequently annihilate in 
the detector, much of the energy of the annihilation products may be detected in the BSC. 
The scatterplot is consistent with these expectations, indicating the two tracks are really 
the proton and anti-proton. 

Fig. 5 shows the distribution of secondary vertices in the xy plane of 7AA candidates in 
the Xcj mass region (error bars). This distribution shows good agreement with the secondary 
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FIG. 3: Mass distribution of AA candidates. Histogram with error bars is data, and the shaded 
histogram is from A side bands events (normahzed). 
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FIG. 4: Energy deposited in the BSC of the proton and antiproton tracks for selected Xcj — > AA 
events. 
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vertex distribution of selected ip{2S) —>■ AA events (histogram), but is significantly different 
from the vertex distribution of ip{2S) ix^ix^ Jj-^, J/ip ^ pp events (stars), where no 
secondary vertex is expected. This indicates the events in the XcJ mass region are real AA. 
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FIG. 5: Secondary vertex distributions. Dots with error bars are for events with their mass in 
the XcJ mass region, the histogram is for selected '<p{2S) AA events, and the asterisks are for 
selected ip{2S) tt^tt^ J/il^, J /ip — > pp events. The dots and histogram are normalized for greater 
than 1 cm, and the normalization for the asterisks is arbitrary. 



A. Remaining backgrounds 

Background from non AA events is estimated from the A mass sidebands as shown in 
Fig. 3, and this can be described in fitting the AA mass spectrum by a linear background. 
The background from channels with AA production, including 'ip{2S) —>■ AA, 'ip{2S) S'^S°, 
^-(25) ^ ATP + c.c, ^{28) ^ E'TP + c.c, V^(2S) ^ iXcJ,XcJ ^ S^SO ^ 77AA, and 
ip{2S) —>■ n^n^ J/tp -n-^n^pp, are simulated by Monte Carlo. By using the branching 
ratios of ^^{28) AA, 'ip{2S) S°S°, and ip{2S) — * t^^t:^ J/ip — > tt^tt~pp measured by 
previous experiments [4], and naively assuming ip{2S) AS° + c.c. and 'ip{2S) — > S°S° + c.c. 
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are one order of magnitude smaller than ip{2S) — > AA and ■0(25') — > E°E°, and XcJ — E°E° 
is about the same as XcJ AA, we obtain the expected total background plotted in Fig. 6. 
The curve in this plot indicates the best fit of the background mass spectrum from 3.2 to 
3.65 GeV/c^. The background from events with more photons is smaller, and Monte Carlo 
simulation of ip{2S) 5°5° indicates that its contamination to the Xc signal is negligible. 
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FIG. 6: Invariant mass distribution of AA selected from Monte Carlo simulated background events 
normalized to the total number of ip{2S) events in the data sample. The curve shows the best fit 
of the mass spectrum below 3.65 GeV/c^. 

B. Fit to the mass spectrum 

Fixing the mass resolutions at their Monte Carlo predicted values ((12.7 ± 0.9) MeV/c^, 
(9.4 ± 0.3) MeV/c^ and (9.8 ± 0.4) MeV/c^ for 

Xco, Xci and Xc2, respectively), and fixing 
the widths of the three XcJ states to their world average values [4], the mass spectrum 
was fit with three Breit-Wigner functions folded with Gaussian resolutions and background, 
including a linear term representing the non AA background and a component described 
in the previous subsection representing the AA background with the global normalization 
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factor floating to take into account possible systematic bias in the background estimation 
(mainly branching ratio uncertainties) . The unbinned maximum likelihood method was used 
to fit the events with AA mass between 3.22 and 3.64 GeV/c^, and a likelihood probability 

of 27% was obtained, indicating a reliable fit. The number of events with errors determined 

I /I O I Q ^ I Q 

from the fit are 15.2_4;q, 9.0_3;"]^, and 8.3_3;4 for Xco, Xa and Xc2, respectively. The statistical 
significances of the three states are 4.5cr, 3.5a and 2.6a. Fig. 7 shows the fit result, and the 
fitted masses are (3425.6±6.3)MeV/c2, (3508.5±3.9) MeV/c^ and (3560.3±4.6) MeV/c^ for 
XcO) Xci and Xc2, respectively, in agreement with the world average values [4]. The detection 
efficiencies from the Monte Carlo simulation were determined to be = (6.07 ± 0.24)%, 
^Xci ~ (6-65 ± 0.25)% and = (6.09 ± 0.24)%, where the errors come from the limited 
statistics of the Monte Carlo samples. 

50 p I 
40 - 




AA mass (GeV/c^) 

FIG. 7: Mass distribution of 7AA candidates fitted with three resolution smeared Breit-Wigner 
functions and background, as described in the text. 
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IV. NUMBER OF iJ^{2S) EVENTS 



The number of ip{2S) events is determined using ip{2S) tt+tt J/ip, J/ip — > pp. There 
are many advantages in using this channel to determine the number of events: 

1. It has the same kind of charged tracks as the channel of interest, and the momenta 
in these two channels are similar, so that in the branching ratio measurement, the 
systematic bias in tracking, kinematic fit, triggering, particle ID, geometric acceptance 
of charged tracks, etc. will cancel out. 

2. It is easy to select, and the error on the branching ratio is small ((2.12 ± 0.10) x 10~^ 
for the world average) [4] . 

The selection criteria of this channel are the same as for the Xcj analysis, except 

the photon is not considered. The invariant masses of 7r~p and Ti^p are required to not be in 
the A mass region to remove '4^{2S) AA background. Fig. 8 shows the pp invariant mass 
distributions of both data and Monte Carlo. There is a huge J/ip signal on top of very low 
background. 

The number of J/ip — > pp events is estimated by subtracting sideband events for pp 
invariant mass regions from 3.0 to 3.05 GeV/c^ and from 3.15 to 3.2 GeV/c^ from the signal 
region {pp invariant mass from 3.05 to 3.15 GeV/c^), giving 

= 1826 ± 44. 

Using the same method, the efficiency is determined using Monte Carlo data as 

e = (17.88 ±0.12)%. 

Using the BES branching ratio for ij;{2S) n+n-J/ip ((32.3 ± 1.4)% [9]) and the PDG 
branching ratio for J/j/j — > pp ((2.12 ± 0.10) x 10"^ [4]), the number of ■0(2-5') events is 
obtained 



"^^^^^ B{ilj{2S) TT+TT- J/^)-B(J/^ ^ pp) 

= (14.91 ± 0.36 ± 1.13) X 10^ 

where the first error is statistical and the second is systematic, including the statistical error 
of the efficiency, the errors from the two branching ratios used, and the uncertainty due to 
the Monte Carlo simulation of the angular distributions. 
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FIG. 8: Distribution of pp invariant mass of tp{2S) —> tt+tt J/V', J/ip pp data (top) and Monte 
Carlo (bottom). 



It should be noted that the efficiency correction factors due to the differences between 
data and Monte Carlo data in the particle ID, the kinematic fit, tracking, etc. are not 
considered, because the same differences exist in the XcJ analysis and will cancel in 

the XcJ A A branching ratio measurement. 

As a consistency check, one can apply the particle ID correction factor (1.043±0.011) and 
kinematic fitting correction factor (0.943 ±0.010), which are measured in following sections. 
One then obtains N^(2S) — (15.16 ± 0.37 ± 1.16) x 10^, which agrees with the number of 
ip{2S) events determined using either inclusive ip{2S) — > Tr'^'!r~J/ip or inclusive hadrons. 
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V. EFFICIENCY CORRECTION AND SYSTEMATIC ERRORS 



The systematic errors in the branching ratio measurements come from the efficiencies of 
the photon ID, particle ID, kinematic fitting, low energy photon detection, MDC tracking, 
the branching ratios used, the number of 1/^(23) events, the A mass cut, etc. 

A. Photon ID 

The fake photon multiplicity distributions in both data and Monte Carlo simulation are 
checked with ip{2S) Tr~^Tr~J/ip, J/ip ^ pp events. The Monte Carlo predicts too many 
fake photons at very low energy (less than 50 MeV) . Using a photon energy cut at 50 MeV or 
reweighting the Monte Carlo events with the measured fake photon multiplicity distribution 
indicates that the Monte Carlo simulates the data with a precision of 4%. This will be taken 
as the systematic error on the photon ID. 

B. Particle ID 

Samples of tt"*", tt", p, and p tracks are selected in ip{2S) tt+tt" J/-^, J/ip pp events 
by requiring a good kinematic fit to this process and good particle identification of the other 
three charged tracks involved. This allows a measurement of the particle ID efficiency, and 
a correction factor of 1.043 ± 0.011 to the Monte Carlo efficiency is found for the channels 
that wc are studying. The error is from the limited statistics of the samples used and is 
taken as the systematic error of the particle ID. 

C. Kinematic fit 

The bias due to the kinematic fitting is caused by differences between data and Monte 
Carlo data in the fitted momentum and error matrix of the charged track and differences in 

the measurement of the energy and the direction of the neutral track and their uncertainties. 
The effect is studied for charged tracks and neutral tracks separately. 
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1. Charged tracks 



The bias from the kinematic fit of the charged tracks was checked using i>{2S) — > 
TT'^TT" J/ip, J/ip ^ pp events. This channel is very clean and can be selected without the 
help of a kinematic fit. By comparing the number of events with and without a kinematic 
fit, the efficiencies for prob^2 > 1% are measured to be (85.14 ± 0.92)% and (90.32 ± 0.24)% 
for data and Monte Carlo, respectively. This results in a correction factor for the Monte 
Carlo efficiency of 0.943 ± 0.010 for this specific channel. 

2. Neutral tracks 

The effect of neutral track measurement is studied using ip{2S) 7XcJ, XcJ 7r+7r~pp 
events. A careful calibration of the neutral cluster information in the BSC (including the 
energy and direction measurement and their errors) was performed using radiative Bhabha 
events from the same ip{2S) data set. By applying this calibration to both data and Monte 
Carlo, the relative changes in the branching ratios of XcJ ~^ 7r~^7r~pp are measured to be 
1.1%, 1.9% and 4.2% for Xco, Xa and Xc2, respectively. No corrections to the efficiencies are 
made; the largest difference (4.2%) is taken as the systematic error in the measurement of 
neutral tracks. 

D. Photon detection efficiency 

The low energy photon detection efficiency is studied with ip{2S) — >• ti~ J / ip , J j'^ ^ 
'K^'K~'K^ events produced in the same data sample used for the XcJ analysis. We assume 
the lower momentum positive and negative charged tracks are the tt"^ and 7r~ from ip{2S) 
decays, and the largest energy neutral cluster is a photon from the 7r° decay. Assuming the 
second photon from the 7r° decays is missing, we do a two constraint kinematic fit requiring 
all the final particles come from ip{2S) decays and the two photons form a n^K The fitted 
four-momentum of the second photon is taken as a test beam into the detector and used 
to determine the detection efficiency. A total of 2901 photons are selected for the efficiency 
study. The same analysis is performed with Monte Carlo events, and agreement between 
data and Monte Carlo data is observed at a precision of 8% for the photons accompanying 
Xco, Xci and Xc2- 
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For converted photons, no specific study was performed since this occurs for only a very 
small fraction of the events (less than 1%), and the difference between data and Monte Carlo 
simulation should be even smaller and negligible compared to the quoted systematic error 
for the photon efficiencies. 

E. Other systematic errors 

The angular distributions of the photon accompanying the Xcjs and the angular dis- 
tributions of the A or A decays may cause a systematic error at the 10% level. This is 
determined by comparing different theoretical models for the angular distributions. The 
uncertainty in the angular distribution of the proton in J/ip decays results in a 4% error in 
the determination of the number of ip{2S) events. 

The Monte Carlo simulated mass resolution may have a bias at the 10% level. This is 
determined from the comparison of A and J/tp signals in various channels involved in this 
analysis. Changing the mass resolutions used in fitting the XcJ mass plot produces small 
changes in the number of events; the maximum change in the three cases is around 3%. This 
is taken as the systematic error due to the mass resolution uncertainty. 

The background estimation, including the uncertainties in the branching ratios used, the 
uncertainties in the simulation of the contamination probability, the parameterization of the 
background shape, and the fitting range used, etc., causes an uncertainty at the 10% level. 
The systematic errors on the branching ratios used, hke B{'^{2S) — > tt+tt" J/^'), B{J/'^ — > 
pp), B{'^{2S) — > 7Xcj) and B(A — > 7r~p) are obtained from other experiments [4, 9]. 

F. Total systematic error 

Table. I lists the systematic errors from all sources, as well as the correction factors to 
the Monte Carlo efficiency for particle ID and the kinematic fitting of charged tracks. Since 
these two correction factors cancel out in the calculation of branching ratios, there are no 
corrections to the efficiencies determined by Monte Carlo simulation for the XcOi Xci and Xc2 
branching ratios, and their errors are not considered in the summation. 
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TABLE I: Summary of systematic errors and the efficiency correction factors. Efficiency correction 
factors are only determined for the particle ID and the kinematic fitting of charged tracks. Since 
these correction factors cancel in the branching ratio calculation, they are not used. 



Source 


XcO Xcl Xc2 


MC statistics 


4.0% 3.8% 4.0% 


Fake photon 


4% 


Particle ID 


1.043±0.011 


4C-fit (chrg) 


0.943±0.010 


4C-fit (neut) 


4.2% 


Phot. eff. 


8% 


Gamma conversion 


<1% 


Angular distr. 


10% 


Mass resolution 


3% 


Background 


10% 


1^(28) number 


8.0% 


B{A TT-p) 


1.6% 


B{i^{2S) ^ 7Xc,/) 


9.2% 8.3% 8.8% 


Total systematic error 


22% 21% 22% 



VI. RESULTS AND DISCUSSIONS 

The branching ratios of XcJ ~^ can be calculated with 



B{xcj ^ AA) = 



Using numbers from above, one gets 

B(xco ^ AA) = {A.7t\l ± 1-0) X 10-^ 

B(Xci ^ AA) = (2.6lJ:° ± 0.6) X 10-^ 
B{xc2 ^ AA) = (3.31} j ± 0.7) x 10-^ 
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TABLE II: Summary of numbers used in the branching ratio calculation and branching ratio results. 
Rb, defined in the text, is the relative branching ratio of XcO to that of iIj{2S) — > tt+tt" J/'0. 



quantity 


XcO 


Xcl 


Xc2 


^obs 


-■^^•^-4.0 




8.3t|I 


e (%) 


6.07 ±0.24 


6.65 ± 0.25 


6.09 ± 0.24 


iV^{25)(10^) 




14.9 ± 1.2 




B{K TT-p) [4] 




0.639 ± 0.005 




B{^l^{2S) ^ jxcj) (%) [4] 


8.7 ±0.8 


8.4 ±0.7 


6.8 ±0.6 


i3(xcJ^AA)(10-4) 


4.7ti| ± 1-0 


2.6l^:^ ± 0.6 


3.3l:}j±0.7 


obs 
"•tv+tt- J/ip 




1826 ± 44 




^TT+TT-J/lP (%) 




17.88 ±0.12 






2.45l[j:i ± 0-46 l-33l[j:|6 ± 0.25 1.33l[j:i ± 0-25 



where the first errors are statistical and the second are systematic. The numbers used and 
results arc summarized in Table. II. 

Compared with the corresponding branching ratios of XcJ PP [4], the branching ratios 
of Xcl a-nd Xc2 AA agree with the corresponding branching ratios to pp within two sigma. 
This is somewhat in contradiction with the expectations from Ref. [3], although the errors 
are large. 

As for XcO AA, the measured value agrees with the pp measurements from BES and 
E835 [2, 10] within 2 standard deviations. One should also note that there is no prediction 
for Bixco ^ AA). 

What we actually measure in this analysis is the relative branching ratio of Xco ^ AA to 
'il^{2S) — > 7r'^7r~ J/i/j. The relative branching ratio is found with the following formula 

^ B{H2S) ^ 7Xcj) ■ Bjxcj ^ AA) • B{A n-pf 
^ B{i){2S) T:+T:-J/i}) ■ B{Jli) pp) 

These results are also shown in Table II. 
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VII. SUMMARY 



AA events are observed for the first time in XcJ decays using the BESII 15 miUion ■0(2-5') 
event sample, and corresponding branching ratios are determined. The results on Xci ^-nd 
Xc2 decays only agree marginally with model predictions. 
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